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Abstract 

Vibration acceleration levels on large space platforms exceed the requirements of many 
space experiments. The Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) is 
being built by the NASA Marshall Space Flight Center to attenuate these distuibances to 
acceptable levels. G-LIMIT uses Lorentz (voice-coil) magnetic actuators to levitate and isolate 
payloads at the individual experiment/sub-experiment (versus rack) level. Payload acceleration, 
relative position, and relative orientation measurements are fed to a state-space controller. The 
controller, in turn, determines the actuator currents needed for effective experiment isolation. 
This paper presents the development of an algebraic, state-space model of g-LIMIT, in a foim 
suitable for optimal controller design. The equations are first derived using Newton’s Second 
Law directly, then simplified to a linear form for the purpose of contiollei design. 
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Introduction 

Acceleration measurements on the U.S. Space Shuttle and the Russian Mir Space Station 
have shown that on-orbit acceleration environments are noisier than once expected [1]. The 
acceleration environment on the International Space Station (ISS) likewise is not as clean as 
originally anticipated; the ISS will not meet its microgravity requirements without the use of 
isolation systems [1], [2], While the quasi-static acceleration levels due to such factors as 
atmospheric drag, gravity gradient, and spacecraft rotations are on the order of several micio-g, 
the vibration levels above 0.01 Hz exceed 300 micro-g rms, with peaks typically teaching milli-g 
levels [3]. These acceleration levels are sufficient to cause significant disturbances to many 
experiments that have fluid or vapor phases, including a large class of materials science 

experiments [4]. 

The Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) is designed to 
isolate experiments from the high frequency (>0.01 Hz) vibrations on the Space Shuttle and the 
ISS, while passing the quasi-static (<0.01 Hz) accelerations to the experiment [51. The 
acceleration-attenuation capability of g-LIMIT is limited primarily by two factors: (1) the 
character of the umbilical required between the g-LIMIT base (stator) and the g-LIMIT 
experiment platform (flotor), and (2) the allowed stator-to-flotor rattlespace. A primary goal in g- 
LIMIT design was to isolate at the individual experiment, rather than entire rack, level; ideally g- 
LIMIT isolates only the sensitive elements of an experiment. This typically results in a stator-to- 
flotor umbilical that can be greatly reduced in size and in the services it must provide. In the 
current design, g-LIMIT employs three umbilicals to provide experiments with power, and with 


data-acquisition and control services [6]. 
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In order to design controllers for g-LIMIT it was necessary to develop an appropriate 
dynamic model of the system. The present paper presents an algebraic, state-space model of g- 
LIMIT, in a form appropriate for optimal controller design. 

Problem Statement 

The dynamic modeling and microgravity vibration isolation of a tethered, one-dimensional 
experiment platform has been studied extensively by Hampton, et al. [5, 7, 8]. It has been found 
that optimal control techniques can be effectively employed using a state-space system model, 
with relative-position, relative-velocity, and acceleration states. In these studies the experiment 
platform was assumed to be subject to Lorentz (voice-coil) electromagnetic actuation, and to 
indirect (umbilical-induced) and direct translational disturbances. 

The task of the research presented below was to develop a corresponding state-space 
model for g-LIMIT. Translational and rotational relative-position and relative-velocity, and 
translational acceleration states, were to be included in the system model. The g-LIMll dynamic 
model must incorporate indirect and direct translational and rotational disturbances. 


System Model 
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A schematic of g-LIMIT is depicted in Figure I. The stator, fixed in reference 
frame 0 , is rigidly mounted to the ISS. The flotor, frame 0 - is magnetically levitated above the 

stator by six Lorentz actuators (two shown), each consisting of a flat racetrack-shaped electrical 
coil, with an active linear (straight) region positioned between a set of Nd-Fe-Bo supermagnets. 
For more information on the basic actuator design see Reference [9], The coils and the 
supermagnets are fixed to the stator and flotor, respectively. Control currents passing through the 
coils interact with their respective supermagnet flux fields to produce control forces used for 
flotor isolation and disturbance attenuation. 

The flotor has mass center F* and a dextral coordinate system with unit vectors / ( , 

and f and origin F 0 . The stator (actually, stator-plus-ISS) has mass center S * and a dextral 

coordinate system with unit vectors s 2 , and X, and origin S 0 . The inertial reference 

frame 0 is similarly defined by n,, n 2 , and tX, and origin N 0 . The umbilical is attached to the 
stator at S u , and to the flotor at F u . When the flotor is centered in its rattlespace (the “home 
position), F* and F u are located at stator-fixed points F,* , and F uh , respectively. 
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Figure 1. Schematic of g-LIMIT ' 


State Equations of Motion 


Preliminaries 


Let E be some flotor-fixed point of interest for which the acceleration is to be determined. 


If E has inertia] position r ? then its inertial velocity and acceleration are i_ N ^ F ^ \-A f 0 £/ 


■d r N d 


r — r 

-No E dt dt \-NoE 


, respectively. (The pre-superscript indicates the reference frame assumed 


fixed, for purpose of differentiation. The first and second subscripts indicate the vector oiigin and 
terminus, respectively.) The angular velocity and angular acceleration of the flotor with respect to 


the inertial frame are represented by ‘^(ti F and s a F , respectively, where a ^ 


(V). 
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fIf* 

Let F be the resultant of all external forces acting on the flotor; M (or simply M), the 
moment resultant of these forces about F ; m, the flotor mass; and £ j (or / ), the cential inertia 
dyadic of the flotor for f , f , and f . Then Newton’s Second Law for the flotor can be 


expressed in the following two forms: 



F = mr No F* 

(Eq. 1) 

and 

M = £• V* V‘x(/ V). 

(Eq. 2) 

From Equation (2), 




N a F = r' \m- Vx(/- v)] . (Eq.3) 


Equation (1) presents the translational equations of motion for the system (in vector form) 
in terms of the acceleration of flotor mass center F* . However, the primary objective of the 
control system design is to attenuate the absolute acceleration of point E, not of point F . 
Consequently, it will be useful to find an expression for r N{)F * in terms of r NoE . The three s t 

measure numbers of r N , will be used as states in the state equations ultimately to be used for 


controller design. 
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In addition to the above objective, the controller must prevent the flotor from exceeding 
its rattlespace limits, i.e., from bumping into the stator. It is desirable, then, that relative-position 
and relative-velocity information (both translational and rotational) also be included in the 
equations of motion for controller design. The s, measure numbers of the relative position and 

relative velocity forF\ from its home position F h , will provide six translational states in the 
system state equations. Six corresponding rotational states will be included as well. The 
remainder of this paper will be concerned with expressing the system equations of motion 
[Equations (1) and (2)] in terms of the selected fifteen states. 


Translational Equations of Motion 

To develop the translational equations, in terms of the chosen states, begin with the 

+ frr* - r v „ c - • (Eq- 4 ) 


following: 


} -F l h r* - Ln 0 e t Lef’ - -N 0 s 0 Ls 0 F h 


Differentiation of Equation (4) yields 


t F,* F* ~ £-NqE — X -EF* 


^NqS q 5? x -.v„f a * 


(Eq. 5) 


A second differentiation gives 


Z-FhF* ~Fn 0 e 


N F , K 

+ a xr EF * + 


- V x 


J F i N F \ 

CO X ^ (0 x Tef* }l 

k s x Vt)' 


r N 0 S 0 ff X Ts () A 


(Eq. 6) 


Substitution for N a f from Equation (3) into Equation (6) yields 


!~F,, F 


= r 


n^e ' 
% S X 


{/"' [m_- % f X (r % F )]} xr tf ‘ - 
( x Ls 0 Fh ) + 1*1 F x ( ® x Lef" )• 


Ln 0 s 0 


- V X r_ s f. 


0 r h 


(Eq. 7) 


In these equations, from the addition theorem for angular velocities, 
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V = V + s co F . 


(Eq. 8) 


Under the assumptions that N co s and N a s are negligibly small and, therefore, that 


Ln 0 s 0 ~ —N a S tl ’ 


(Eq. 9) 


Equation (7) reduces to 

z!f a v = ryv 0 £+U 1 


M- S (ti r x 


(/•^w F )]} ~LnqS ii + 5 w f x( S w F xr £f *). (Eq. 10) 


Linearization about ftL =0 yields the following result: 


>!/.,>* = iiw 0 F + k"' 'Ml 

NS 


x -i:r* Ln q s u 


(Eq. 11) 

(Note that assuming 'co J to be negligible does not imply that reference frames (0 nnd © aie 


identical; it means rather that (s) can be treated as if it is in pure translation relative to foi 
the frequencies of interest.) 

Equation (11) describes in terms of the acceleration of an arbitrary flotoi -fixed point L. 


For E located at F*, Equation (11) can be used straightforwardly with Equation (1) to yield 


rn 






\E In,a 


O^w 


(Eq. 12) 


where the second term on the right-hand-side is the indirect disturbance acceleration. Define now 
an (unknown) indirect translational acceleration disturbance input to the flotor, applied at the 
stator end of the umbilical: 

a -r (Eq- 13) 

— in — —NqS u K 4 

Substitution from Equation (13) into Equation (12) yields the following: 


r -V 


m 




E ~ a in . 


(Eq. 14) 
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Turn next to Ln 0 e- Substituting from Equation (14) into Equation (11), using Equation 

(13) , and solving for one obtains the following: 

1 3 , 

Ln E~ — E + L ' ■ Eq. 15) 

m J ' 

Appropriate expressions for F and M will now be determined in terms of the position and 
velocity of the flotor, relative to the stator. These expressions will be substituted into Equations 

(14) and (15) to obtain more useful representations of the equations of motion. 

The force resultant F is the vector sum of the six actuator (coil) forces F'. (/' = 1,...,6) , with 

resultant F r ; of the umbilical force F u , caused by umbilical extensions and rotations from the 
umbilical home position; of the (unknown) direct disturbance forces, with resultant F ( /; and of the 
gravitational force F s . Gravity may be neglected for a space vehicle in free-fall orbit. (Gravity 
gradient and fluctuations in the gravitational field are lumped together in F ( /). The moment 
resultant M is the vector sum of the various moments acting on the flotor. These moments are 
from four sources: the coil forces, the umbilical forces, the umbilical moment, and the direct 

disturbance moment. Let M r represent the moment about F , due to the respective coil forces 

£ 

F t , assumed to act at their respective coil centers B r Represent by M_ uf the moment about F 
due to the umbilical force jF u . Let M_ u represent the umbilical-transmitted moment, due (1) to 
translations of umbilical attachment point F u from its home position F uh , (2) to rotations about 
F u , of the flotor relative to the stator, from the home orientation, and (3) to umbilical bias 
moment in the home position. Finally, let M_ (l represent the (unknown) moment due to the direct- 
disturbance forces and moments, with F_ d assumed to act at F* . Tn equation form, 
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£=£,+£. +£j. 

(Eq. 16) 
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where 

34 

w- 

II 

I- 

(Eq. 17) 

and 

= +M uf +M U +M d , 

(Eq. 18) 
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where 

Mr = X — < ’ 

i=l 

(Eq. 19) 

and 

M«/ —Lf*f„ X — h • 

(Eq. 20) 


More explicit expressions for F|. and F (( will now be developed. If the actuator has coil 
current length L, , and magnetic flux density B l B t , then the associated actuator force 


becomes 


£' =-I l L i B i l i x| r (Eq.21) 

Assume a translational stiffness K\; j in the 5^ direction for an umbilical elongation in 
the s_j direction, and a corresponding translational damping C l ,j J . Assume also a translational 
stiffness K\’J in the Sj direction for a positive umbilical rotation about the 5 . direction, and a 

corresponding translational damping . Let F_ b represent the umbilical bias force, exerted by 
the umbilical on the flotor in the home position. Then the total force of the umbilical on the flotor 
becomes 


3 3 


L u = 


£ IkA 




3 3 


S: \S ; + 


M M 
3 3 


+ 


£ £*i J 

7=1 i = 1 


FfS 


9 * S; 


£ £c' 

j = 1 /=• 

3 3 

lly + £ £c 

7=1 f=l 


dr 

s d_ 

dr 




\*J 


(Eq. 22) 


he-l,)-!, 


& 


r + Kb > 
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where the vector F/S Q represents the rotation about F u of the flotor, relative to the stator, from the 
home orientation. Equations (17), (21), and (22) now provide expressions for F ( . and F„ that 
can be used in Equation (16) for substitution into Equations (14) and (15). 

As with F', and F u above, M_„ can also be expressed in more explicit form, in analogous 

fashion. Assume a rotational stiffness K' r '/ and a rotational damping C r '/ , about the diiection, 
for umbilical twist about the £. direction. Assume also a rotational stiffness K' r ’/ and a rotational 
damping C' r / , about the ^direction, for umbilical extension in the s, direction. Let M_ b represent 

the total umbilical bias moment, exerted about F* , by the umbilical on the flotor, in the home 
position. Then the moment M_ u can be expressed by the following: 



Substituting Equation (21) into (19), and Equation (22) into (20) produces expiessions for M_ ( . 
and M_ uj respectively. These expressions, along with Equations (17), (21), (22), and (23) can be 
substituted into Equation (15) to produce a more useful expression for Ln^e- 

In principle, one could select a set ot umbilical elongations, elongation lates, lotations, and 
rotation rates as the system states. However, since the translational form of Newton’s Second Law 
given by Equation (14) employs the acceleration r F * F * , it is more natural to choose states based 

on r r *r* rather than r F F .To accomplish this one must substitute appropriately for r FF and 

— t h f — r uh r u 

its derivative into Equations (22) and (23). The pertinent relationship is 
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r .. r = r * + r m ,+r, 

h " hF “ F uh F h F h F F F u 


(Eq. 24) 


Note that r . and r , are constant vectors of equal magnitude and opposite direction in the 

~ F uh F h ~ F F „ 

stator- and flotor-fixed reference frames, respectively; in particular, 


(Eq. 25) 


It will be shown later that Equations (24) and (25) facilitate expressing the umbilical foices and 




-L 

r -F: hF ; -h 

> — « 

- r - F ' Fu 'L f 

' —3 ^ 


~tF*F u -/ 3 


„ ' d ( \ . F/S a S 

moments in terms of the following twelve scalar quantities: f f * f * -s t , )' -i » 2 -< ’ 


anc l ( F/ ' v 0)-|. : . These quantities will be chosen as system states, for the state-space 

dt v ' 

description. It is also desired to include, as auxiliary states, the measure numbers of LnqE [ see 

Eq. (15)], for control purposes. Then Equations (14) and (15) will yield the translational 
equations of motion in state-space form, provided that (as will be shown later) and r 




be expressed in terms of -{r F * F * ) and respectively. To this end, recall that 


can 


dt 


dt 


% s and N a s have been assumed to be negligible in the frequency range of interest; therefore, 
one has the following: 


r * * 

-FhF 




dt 


(Eq. 26) 


hu- =^rh,:r-) 


and 


(Eq. 27) 



Similarly, 


_ S d 2 / \ 

-—T^NnEj 


— NqE ^ 2 i 


(Eq. 28) 


Note that although the terms involving N (ti S and V 0C 5 have been neglected in Equation (27) their 
effects in Equation (14) could alternatively be considered as lumped together with Ln 0 s u ' n 2 m • 
Including these higher-order terms, Equation (13) would be rewritten as follows: 

■^knr'Y V«(Vxr fl y). (Eq. 29) 


2,„ = W,- Vxr f ; f - 2 'to' x 


NS 


Equations (14), (16), (17), (21), (22), (26), and (27) provide the basis for a state-space 
form of the translational equations of motion, using as states the 5. components of Lf^f* ' 




FIS 


0, and — ( r/5 e). In addition, for control purposes, low-pass-filtered 


dt 


approximations of r N[)E [beginning with Eqs. (15) and (28)] will be used as auxiliary states. This 
development will be presented later [see Equations (65) through (67)]. 


Rotational Equations of Motion 

The rotational equations of motion are given by Equation (2) [or, equivalently, by 
Equation (3)], where M is defined by Equations (18), (19), (20), and (23). Consider the left- 

hand-side of Equation (3). It will be shown that N a F = — r 'ej, under the assumption that 

at* 

V ancl N u S are approximately zero. 

Differentiating Equation (8), one can express the angular acceleration *V , in terms of 



14 




dt K ~ ' dt v “ 


(Eq. 30) 


or, expanding, as 


,v a F 


J/ 


v)+ ^-(vl+'Vx V- 

“ ' dt. ~ 


(Eq. 31) 


Assuming, as before, that =0 and oc ~0, Equation (31) simplifies to. 


v-£(v). 


(Eq. 32) 


The rotation vector F,s Q, defined previously, represents the rotation of the flotor, relative 
to the stator, from the home position. (Conceptually, F/5 0 is a free vector <t> > wheie 

F/s indicates the positive direction of the rotation axis; and (J), the rotation angle about that 


axis.) The angular velocity vector S (0 F is related to the rotation vector as follows [10]: 

v = i(w 9 ). 

dt v 

Substitution from Equation (33) into Equation (32) yields the following: 


N a F 


=^C'e). 


Substitution from Equation (34) into Equation (3) yields: 


(Eq. 33) 


(Eq. 34) 


^( w s)-r-k-vx(/-vl 

dt " 


(Eq. 35) 


Assuming, as before, that V = 0 , one obtains: 


J -( / r / 5 e )=/ _1 • 


dt 


M_ -V>4 s a F 


(Eq. 36) 
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Linearization of Equation (36) about 'W' = 0 yields the following form of the rotational 
equations of motion: 


dr 


( F/S e )=/ _1 M_. 


(Eq. 37) 


Substitution from Equation (18) into Equation (37) yields: 

T ( F/ 'e)= r 1 .(w, +m„ +m,„ +M„), 


V 


dt‘ 


(Eq. 38) 
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where (from before) M_ c - ^ r F * R . xF|., 

/=! 


(Eq. 39) 




3 3 


3 3 




;=i /=i 

3 3 


;=i i’=i 

3 3 


+X + S S'V j( r,s <>k 


7=1 <=1 


7=1 '=1 


— 7 




(Eq. 40) 


+ Ki, , 


and 


— «/ — — F*F U X — h ■ 


(Eq. 41) 


Equations (38) through (41), along with (21), (22), and (24) through (26), provide the basis for a 
state-space form of the rotational equations of motion. The states are those defined previously. 


Equations of Motion in State-Space Form 

Define the following relative-position states: 

SO that, -FhF* = X nl ll + X a2$-2 + 'E/3^3 • 


(Eq. 42) 
(Eq. 43) 


Define next the relative-velocity states, x bi = x aj . 


(Eq. 44) 
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Differentiation of Equation (43), along with the use of Equations (44) and (26), leads to the 
following: 

= x bll\ + x b2$-2 + x b3 §. 3 ' ^) 

A second differentiation, along with the use of Equation (27) yields: 

L F*F* = + x b2§-2 +'^Z)3^3- ^4- 46) 

Introduce the use of a pre-superscript in parentheses to indicate the coordinate system used for 
componentiation of vectors. Then Equations (43), (45), and (46) take the following respective 
matric forms: 


(5) 


lf; f * = < 




r = > 


v «3 


(S) 


r F* h F* = —a ~ —b ’ 


and 


(S) 


LfIf' = x b • 


where x a and x b are defined as indicated. 


Define the relative-angular-position states, 


**- w *h. 


so that 


F S 0=^lIl + *dli2+ x <nh ' 


Define next the relative-angular- velocity states, 


x ei - X di • 


d If 


dt 


( F,s e)=^,l 


1 JrX e2$-2 X e 3-3 ’ 


(Eq. 47) 


(Eq. 48) 
(Eq. 49) 


(Eq. 50) 


(Eq. 51) 


(Eq. 52) 


so that 


(Eq. 53) 
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and 


dr 


( f/S 0)= ^,5, +X el $ 2 +-V ( ,3^3 


(Eq. 54) 


The respective matric forms of Equations (51), (53), and (54) are then 


(5) F/S 


0 = 


l d i 


b/2 


= * ,/ . 


1^/3 J 


(Eq. 55) 


(S) s 




(Eq. 56) 


and 


(5) s ;2 , x 

f( F/s s)-,. 

dt~ 


(Eq. 57) 


where x d and x^are defined as indicated. 


Equations (14) and (46) can be used to develop a state-space equation for x b . First, 
express Equation (46) in measure-number form: 


{S) 


r -F>* = < 




l /)2 




: ' 


(Eq. 58) 


Next, represent by the direct translational acceleration disturbance to the flotor, due to 
unknown direct-disturbance force F _ d . In particular, 


m 


a d =-E d . 


(Eq. 59) 


Then Equation (14) can be expressed as follows: 


3V 


_1_ 

m 


{S) F r + (S) F u 




(Eq. 60) 


Equations (38) and (54) can be used to develop a state-space equation for x e . First 


express Equation (54) in measure number form: 



18 




(Eq.6l) 


Next, represent by a d the direct rotational acceleration disturbance to the flotor, due to unknown 
direct-disturbance moment M_ d . In particular, 


&,= L 'ZLr 


(Eq. 62) 


Then Equation (38) can be expressed as follows: 


i,= <s '{r' ■(M-, +K,+K,„)+a <i y 


(Eq. 63) 


Define a rotation matrix S * F Q that describes a coordinate transformation from the flotor-fixed to 
the stator-fixed frame, so that 


: A = \ s,f q] £ 


(Eq. 64) 


Then Equation (63) can be re-expressed as 


i= s/F e{' 


1 f, r' [ ,n M I , + ,F) M„ l + "'M_, ]+ " (/,, } 


(Eq. 65) 


( F \ * 

where ' 7 is a central inertia matrix of the flotor, for / , / , and / . 


/ r \ 

To approximate r N t using states, define x r by 


® h -NqE ~ —c ’ 


(Eq. 66) 


for some high value of circular frequency co h . Taking the Laplace transform, 


l°>* j J 


(Eq. 67) 
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so that ^ r = (5) lN 0 E > (Eq. 68) 

for to«o),, . 

Equation (15) can be expressed, in expanded form, as 
Im.£= — {F.,. +F U ) + 1 ' +M_ U + M_ uf )xr F * E +a d xr F » E +a d . (Eq. 69) 

w 

Note that Equation (69) is a vector equation; to express it in measure-number (i.e., matric) form, 
one must employ a matric notation for the cross product. To this end, observe that, for arbitrary 
vectors 

L\ = x \f_\ + >'iZ 2 + z 'Li (Eq. 70) 
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Six of the state equations for the system are given by Equations (48) and (56); nine more, 
by Equations (60), (65), and (74). The latter nine are written in terms of the various forces and 
moments acting on the system. The force loads are defined in vector form by Equations (17), 
(21), and (22); the moment loads, by Equations (39) through (41). In the following section these 
loads will be written in measure-number form and expressed in terms of the states, for substitution 
into Equations (60), (63), and (74). 

Force Loads 

In matric form, the i"' control force (Eq. (21)) can be expressed as 

- f- 1, <S t S!F Q B '% ]/, . (Eq. 75) 

Define the control input w, , associated with the i lh control force, by 

(Eq. 76) 

The resultant control force becomes, 

( 5 ) £r=E (S) £r =F c“i (Eq.77) 

/=! 

where F r and u are a row vector and a column vector, respectively, defined as indicated. 

Introduce the following skew-symmetric, small-angle representation for the coordinate 
transformation matrix [10, (p. 352)]: 

■^ Q= (/+W' r/ 'V), (Eq. 78) 

where / is the identity matrix. Then, substituting Equation (78) into Equation (75), one has 

(5> £;. =[- L, iS) ti (/+ (S)r/5 e x )fl, (F) S,l/,. (Eq. 79) 
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Expanding, 


w £' 


<*>/* (F) g _ {S) j* (F) g X (S)F/.?0 


(Eq. 80) 


Next, substitute into Equation (80) from Equation (55), and represent the current I { as the 
sum of a bias value /„ and a fluctuating value I F . This yields the following expression for the 


control force: 


( cx ; (*^) A X ( /"") A / \ (5) A X ( F") ~ X 

<s, £' = -L, B, /, "/, B, + L, B, (/„, + / f J /, B, x, . 


(Eq- 8 1 ) 


The bias current is the current portion necessary to counteract the bias force and moment required 
to hold the flotor at its home position, and the fluctuating current is the additional portion 
necessary for control. Assuming that the fluctuating current and the rotational states are small, 
Equation (81) linearizes to the following: 

( S) ; x ( F )- x ' 


(SUX ( F) a 

Li 


‘V. =\- lb , -t; •' % /,+ / «/;, 7. B, A,. 


(Eq. 82) 


The resultant control force can now be expressed as follows: 


(S) 


E r ='^ ,S) F l ,=F,,,u+F,„x„, 


1=1 


where (as before) 


w, = T i - 


and F ni , and F vd are defined as indicated, from Equation (82). 


(Eq. 83) 
(Eq. 84) 


Turn next to the umbilical force. Equations (22) and (23) can be expanded into matric 
form, in terms of the selected state variables as follows. Recall from Equation (25) that 


(S). 


F tth A 


- (F) r . 

~F F u 


(Eq. 85) 


Then Equation (24) can be represented as follows: 
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(S) 


r F , = <s) r . .+(’' F fi-/) <f, r . . 

- F «/i F i< _ F A F ' y F F„ 


(Eq. 86) 


Assuming small angles Equation (86) can be re-expressed as follows: 


($)- _(S)_ (S) F/S fl x (F) 

f TT TT E / 


-F uh F u 


■F h r 


~F m F„ ’ 


(Eq. 87) 


or, in terms of the states, 


(S) (F) x 

' r, , =.r„ - r , 

-6,/iA - a “F F„ 


(Eq. 88) 


Si i 

One can also represent ~~~\^.f hF J > n terms of the states, as follows. First differentiate Equation 


dt 


(24), with the stator frame assumed fixed: 

f )=- 
(It ^ /,F " ' (It 




r , . 

“ Ft, F 
V h J 


S F 

+ co xr * 
- 'Z 7 


(Eq. 89) 


Substituting from Equation (34) into Equation (89) yields: 




F )=- 

dt ' “ h l< ' dt 




* * 

v F * F , 


+ (W 0 )xr . , 

dt V ~ F F„ 


(Eq. 90) 


or, in terms of stator-fixed coordinates, 


[S)S 


£1 

dt ^ F " hl " 


( s ) s 


dj 

dt 


\ 


- f * p * 

F h F 


+ 


V rhr J 


(5) S 


— ( f / 5 e) 

Jr v 


s,f q tn L 


F F u 


(Eq. 91) 


Using Equation (78), Equation (91) can be re-expressed as 

(S)s At N {S)5 ([ f \ r<«s 


±(. 

dt ^~ F " hF " 


dt[~ F h F 


+ 


^-( //s e)| (/+ (S)F/V 0 x ) (n z; fV • (Eq. 92) 


Substituting with state variables, 


iS)S £ 

dt 


) _ -b + & ( 7 + -<> ) ~F*F U ’ 


(Eq. 93) 


or, assuming the states are small, 

(S) s 


d ( )_ V . (F) _x 

^F uh F u j ±l> 

Finally, the umbilical force, Equation (22), can be expressed in terms of the states by substituting 
from Equations (88) and (94): 

’>£„= i} C„x, + m E>, (Eq.95) 


dt 


r . x 

F F u 1 


(Eq. 94) 


(S) 


where 


and. 




K? 
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K? 

K? 
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1 , 2 
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II 

to 

C 2 ' 1 

r 2,2 

*~7/ 
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^7/ 


Cj‘ 

^3,2 

^ tt 




k ]; 2 

^'; r 

tl 

F 2 ; 1 

t^2,2 

F 2 ; 3 



k?; 2 

F 3 ; 3 


C M 

cjr 

c 1 - 3 ! 

II 

C 2J 

L /; . 

C 2 ’ 3 


C 3J 

K -'tr 

c 3 - 2 

C 3 ’ 3 

v "7r 


( 5 ) 


(Eq. 96) 


(Eq. 97) 


(Eq. 98) 


(Eq. 99) 


Collecting terms, one obtains the following: 

> £,=W„ i „-C„x t -K-Ar„ < V r , k, -{c„-c„ k + m £ t , (Eq. 100) 


or 


(5> £ m = £««*« + + F Md*d + + KJ ’£b . 


(S). 


(Eq. 101) 


where F uta , F H „,, F ((/(/ ,and F Mfe are constant matrices defined as indicated. 
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Moment Loads 


The i' h control force, F'. , exerts on the flotor a control moment, 


Ml. =r F v xf;. 


Substituting from Equation (21), 


M<- ~LF*Bj x ( 1 1 Lj Bj / • x ) ; 


or, in measure-number form, 


(F) M;.=(-/ ( L ( .B, : ) <F) r 


(F) x 


F A 


f,s e <I, i,] <f> s, 


Reversing the order of the second cross product above, and using the small-angle a] 
to v IF Q [Eq. (78)], along with its orthogonality property, one obtains the following: 


cn 






Expanding, and reversing the order of the cross product r/s Q x ( ^ , one has 


(F) -x ( 5 ) a 


i + (W) (n £^ 


(F) a X / ( 5 ) ax 

B, L, x d 


Since 


1 1 I B; ^ F, » 


and since 7 f and x rf( have been assumed small, Equation (106) reduces to 


(f ’£ ISI 1,]+ (/», I,®, <r ’£( 


in (5) 

B \ L x,, 


or 




M = 


LA in r X F % in £ iS) L 


I + 


7 7 77 (H X (H-X(S)aX 

)B, Lf'b, —i Li 


Equation (109) can be written alternatively as 




(Eq. 102) 
(Eq. 103) 

(Eq. 104) 
roximation 

(Eq. 105) 

(Eq. 106) 
(Eq. 107) 

(Eq. 108) 
(Eq. 109) 


(Eq. 110) 
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for M' n and M[ d appropriately defined. Then the resultant moment due to the coil forces can be 


expressed as 


(F> 


o 

M, = Yj M ;. u I, + M\. d x d = M, u u + M rd x d , 


1=1 


where (as before), 


«/ = h » 


(Eq. Ill) 
(Eq. 112) 


and M ai and M rd are defined as indicated. 

The umbilical force F u exerts a moment M_ uj on the flotor, as defined by Equation (20). 


An expression for M_ uf is obtained by appropriately incorporating (S) F„ from Equation (101) 


into Equation (20): 


+ F„„x, ] 


(Eq. 113) 


The moment r x F , , due to the bias force F . , has been omitted from Equation (113) since it 
-Ff„ - b 

is considered part of the bias moment M_ b in Equation (23). Upon linearizing fL Q as before, 
one can rewrite Equation (1 13) as 


( F )*m (F) X 

M r 


Zi-vT 


■ FT, 


(/- (S) F/S e x ) [F ma x tt +F Mb x b + F Md x d +F ute X'] . (Eq. 1 14) 


Equation (1 14) can be simplified under the assumption that the states are small quantities: 


(F) 


M m/ = <F) 4v \ F mX.a+ F «bXb+ F «dZd ]• 


(Eq. 115) 


Substitution of the definitions for F uta , F utb , F utd , and F uh , , given by Equations (100) and (101), 


into Equation (1 15), yields 


% = «v 


~1’F„ 


-K„x„-\ K„-K„ {F V rF \x„-C„x b -\ C„-C„ \x. 


, (Eq. 116) 
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or, expanding. 


(F) 


M »f = 


- (F V. 

“F*F„ " 


( T ) X in 

- r C, ( 
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_(F) r x 


" FF, 


/T, - 
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(H„x 

F*F„ 


_(F) r x 


“F*F„ 


C,.- <n r\ C„ <V 


7T, 


-F*F„ 






(Eq. 117) 


The vector M„ can also be expanded into matric form as a function of the state variables 
by substituting into Equation (23) from Equations (88) and (94), to yield 

( \ ( 


(F W. = ™g 


v rr±d 


where 


and 


Use of Equation (78) with Equation (118) yields 
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, (Eq. 118) 


(Eq. 119) 


(Eq. 120) 


(Eq. 121) 


(Eq. 122) 


. (Eq. 123) 


Assuming that the each of the state variables is small, 


Combining Equations (1 17) and (124), one obtains 


(Djl# 




(F)x 


K„ |x.. + | -C rt - (F V C tt lx,. 


+ I K n {F Y , - - (F V X . 7f (r + (F V X . K„ (F V X „ + (S) M x 1 x, 

" ~FF„ ,r ~r r„ ,r ~F F„ " ~F F„ — 


(B). 125) 


+ 1 c„ tF Y. 

" ~F F„ 


C - (F V X C + ir V C ' r 


(F) x 


(F) x 


-F*F, 


"// -IT* 


FT,, 


1, 




Or, 

(F) M „ + (F) M uf = M uM x (l + M utb x b + M ur(/ x (/ + M ure x t + iS) M b , (Eq. 126) 

for M um ,M llth , M ltrd and M un , appropriately defined. 


State Space Equations 

The force loads ( F_ c and F _ u ) and the moment loads (M_ f . , M _ u , and M _ uf ) have now been 

expressed in measure number form [see Eqs. (83), (101), (1 1 1), and (126)] in terms of the state 
variables. These expressions can now be used to obtain the final, expanded form of the linearized, 
state-space equations of motion. 

Beginning with the state equation for x b , substitution from Equations (83) and (101) into 
Equation (60) yields 


27 = 


1 u (a 

\ m j 


| *a + 


n „ ^ 
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mb 
K m J 


!<> + 


—{Kki+Kd) 
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1 

(±F 1 

1 ute 

!<■ + 

(If ) 

J 

{”> J 


l m J 


u- ( % n +^Q d +-E b -( Eq. 127) 
m 


Recall that the control, u, is a column vector of actuator currents /., each of which is the sum of a 


bias value I B and a fluctuating value l Fj . Then one can define a new control u such that, 
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u=[ B +u, (Eq. 128) 

where [ B is a column vector of actuator bias currents and u is a column vector of fluctuating 
currents. Since the bias current is defined as the current necessary to counteract the umbilical bias 
loads, substitution from Equation (128) into Equation (127) yields the following linear state 
equation for x b : 
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ri F i 
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jn 
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(Eq. 129) 


Turning now to the state equation for x e , substitution of Equations (111) and (126) into 
Equation (65) results in the following expression for x e : 

x= s/F Q (n I ' M ma x a +M utb x b +(M urd +M ( . d )x d +M ure x e + (S) M b +M ru ii + S/F Q iF) a (t , (Eq. 130) 
or, in terms of the new control u * , 

x= s ' F Q ,r) r' + M ulb x b + (, M„ rJ + M„, )x d + M„,x, + M„,u } + S,F Q ,n a, . (Eq. 1 3 1 ) 
Using Equation (78), and combining terms, one obtains 


£= (/+ +M„ b x b + {M, rb + Mjx„ + M ur ,x, +M„ u ] 


+(l+ (s)F,s f ) (n a, . 


(Eq. 132) 


Assuming small state variables, small control, and small direct angular acceleration disturbances, 
Equation (132) reduces to the following linear state equation for r c , : 

Turning now to the state equation for x c , substitution from Equations (83), (101), (111) 
and (126) into Equation (74) results in the following expression for x e : 
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Using Equation (78), one obtains 
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(Eq. 134) 


' (Eq. 135) 


‘j ' (Eq. 136) 


Assuming small states, small direct disturbances, and small controls, Equation (136) reduces to the 
following linear state equation for x r : 
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A state-space representation of the system is given by Equations (48), (129), (137), (56), 
and (133), for state vector 

* = [*! d d d if. (Eq 138) 

and control vector 

IL = |/f, I r 2 I f 3 h A I f 5 h (l ] • •^ t l- ^39) 

In summary, the final form of the linearized state-space equations of motion are as follows: 
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(Eq. 141) 
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(Eq. 142) 
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(Eq. 143) 


kIsA^mJ x b + [ (F) r'{M ur<J +Mj} *A (F) r l Mjx e +[ (n r'M n }u + ( % . (Eq. 144) 
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Verification of State-Space Equations 

The linearized state-space equations of motion for g-LIMIT, given by Equations (140) 
through (144), were verified by comparing the linearized-system open-loop response to the open- 
loop response of a nonlinear “truth” model, for various inputs. Two input scenarios were used for 
these comparisons. In the first scenario (Scenario # 1), the flotor was excited via an actuator by 
an arbitrary force-plus-torque pulse combination, to result in a fully coupled six-degree-of- 
freedom response. In the second scenario, the flotor was excited via an actuator by a force pulse 
in the +Z stator-fixed direction. For each scenario the motion of the linearized system was then 
compared with the response of the nonlinear model. It was found that the nonlinear and linearized 
models exhibited essentially identical responses. For each scenario, the potential and kinetic 
energies were plotted as functions of time. When the systems were made conservative (i.e., with 
the umbilical damping removed) it was shown that the total energy was constant, as expected. 

The nonlinear model was constructed from the nonlinear form of the equations of motion, 

Equations (1) and (3). Employing the quaternion representation (] of the transformation from 

stator-fixed to flotor-fixed coordinates, the following kinematic equation was used propagate the 
rotational motion of the flotor [12]: 


where 
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(Eq. 145) 


(Eq. 146) 


and 


(Eq. 147) 


Equations (16) through (23) were used, in their respective nonlinear vector forms, to compute the 
forces and moments in the nonlinear model. The small-angle approximation 

h /5 0- 5, = 2q,q 4 (( = 1,2,3) (Eq. 1 47 a) 

was used in the third summation terms of Equations (22) and (23), where the factors q, are as 

defined by F S( [ = \ ( 1\ Qi % • (Eq. 147b) 

(This is a reasonable approximation for the expected range of angular motions of the flotor relative 
to the stator.) 

Representative flotor and umbilical parameters, shown in Table l, were used foi the 
simulations. There are three umbilicals included in this model of g-LIMTT. The translational and 
rotational stiffness matrices for each umbilical were assumed to be diagonal along an umbilical- 
fixed set of coordinate directions. Similarity transformations of these diagonal matrices weie 
performed assuming a coordinate transformation from each local umbilical-fixed reference frame 
to the stator-fixed frame. First, a coordinate rotation about the stator-fixed +Z axis of 120 deg and 
240 deg was performed to align umbilicals #2 and #3, in their respective home locations. Then, for 
each umbilical, a 20-deg rotation about each coordinate axis was used to represent an arbitrary 
misalignment of the diagonal-stiffness directions to the stator-fixed directions. The translational 
and rotational damping matrices were assumed to be proportional to the stiffness matrices with a 
damping ratio of 3% used for all of the vibrational modes. The resulting umbilical stiffness and 
damping matricies, given below by Equations (148) through (159), were used in the simulation 
study. The superscript-in-parentheses notation denotes the umbilical identification number for 
each of the three g-LrMIT umbilicals. All stiffness and damping translation/rotational cross-terms, 

i.e., KJ" , K r f'\ C„ (0 , andC rf (0 , were considered to be zero. 
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28.49 -2.46 8.30 

K" } = -2.46 26.73 -5.85 N/m, 

8.30 -5.85 44.77 

' 25.04 0.47 0.91 " 

Kj 2) = 0.47 30.18 10.12 N/m, 

0.91 10.12 44.77_ 

" 29.30 1.99 -9.22" 

K™= 1.99 25.92 -4.27 N /m, 

-9.22 -4.27 44.77 

1.26 -0.04 0.16 " 

C„ (I) = -0.04 1.22 -0.11 N/ m/sec, 

0.16 -0.11 1.58 

1.19 0.00 0.01 ' 

C„ UJ = 0.00 1 .29 0.20 n / m/sec, 

0.01 0.20 1.58 

1.27 0.03 -0.18" 

Cf= 0.03 1.21 -0.08 N/ m/sec, 

-0.18 -0.08 1.58 

1.86 0.09 -0.33" 

K rr m = 0.09 1.93 0.23 N-m/rad, 

-0.33 0.23 1.20 

1.99 -0.01 -0.03" 

Kj 2) = -0.01 1.79 -0.40 N-m/rad, 

-0.03 -0.40 1.20 


(Eq. 148) 


(Eq. 149) 


(Eq. 150) 


(Eq. 151) 


(Eq. 152) 


(Eq. 153) 


(Eq. 154) 


(Eq. 155) 
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and 



1.82 -0.07 0.36 
-0.07 1.96 0.17 
0.36 0.17 1.20 


N-m / rad, 



0.055 0.001 -0.012 

0.001 0.063 0.010 

-0.012 0.010 0.033 


N-m / rad/sec, 



0.062 0.002 -0.002 
0.002 0.056 -0.015 

-0.002 -0.015 0.033 


N-m / rad/sec, 



0.059 -0.004 0.014 
-0.004 0.059 0.005 

0.014 0.005 0.033 


N / m/sec. 


(Eq. 156) 


(Eq. 157) 


(Eq. 158) 


(Eq. 159) 


In addition to the parameters listed in Table #1 the actuator currents were set to initial bias 
values. These bias currents were required to produce a bias force and moment to move the flotor 
from its assumed relaxed position to the home location. The flotor relaxed-position was assumed 
to be 2 mm from the home-position, and to be misaligned by approximately 2 deg about each 
stator-fixed coordinate axis. This resulted in the following set of bias current values: 


7^= -0.264 A, 7fl 4 = -0.1593 A, and 0.1 23 A. 


In Scenario #1 the six actuator current pulses were chosen to approximate combined 
impulsive force loads in the +Z stator-fixed direction, and moment loads about the +X and +Y 
stator-fixed directions. These loads were sized to produce a total energy output of 0.003 N-m 
equally divided among the three loads. The current pulses were shaped as one-minus-cosine 
functions to provide easily integrable inputs with no temporal discontinuities. The pulses were 
initiated at 1 sec and lasted for a 0.1 -sec duration. Table #2 shows the integrated values of the 
actuator current pulses for Scenario #1. 

Figures 2 through 6 show comparisons of the simulation results for the linear state-space 
model and the nonlinear “truth” model for Scenario #1. As shown in these figures the linear model 
responses match the truth model very well for all of the state variables. The state responses evolve 
in time with the expected frequency and damping. The errors between the linear responses and the 
truth-model responses are shown in Figures 7 and 8. These relatively small errors, assumed to be a 
result of the linearization process, diminish as the system response decays. Figure 9 shows a plot 
of the flotor kinetic and potential energy for Scenario #1. As expected, the impulsive loading 
results in an initial translational kinetic energy of 0.001 joules and a rotational energy of 0.002 
joules. The kinetic and potential energies oscillate out-of-phase, and the total energy dampens out 
exponentially. 

A second scenario (Scenario #2) was simulated to demonstrate that the model conserves 
energy when the damping matrices are all set equal to zero. In Scenario #2 the six actuator-current 
pulses were chosen to approximate an impulsive force load in the +Z stator-fixed direction only. 
The impulse moments used in Scenario #1 were eliminated from the actuator loading. As in 
Scenario #1 the loads were sized to produce a total energy output of 0.003 N-m. In this case the 
total energy input to the flotor was along the stator-fixed +Z direction. The current pulses were 
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shaped as one-minus-cosine functions to provide easily integrable inputs with no temporal 
discontinuities. The pulses were initiated at 1 sec and lasted for a 0.1-sec duration. Table #2 shows 
the integrated values of the actuator current pulses for Scenario #2. 

Figures 10 and 1 1 show a comparison of the flotor output energy responses for Scenarios 
#1 and #2, respectively. The potential and kinetic energy responses differ between the two 
scenario, due to the difference in the impulsive loading. However, the total flotor energy remains 
constant at the expected value of 0.003 joules for each of the scenarios, after application of its 
respective impulsive loading, at 1 sec. This demonstrates that the model conserves total energy at 
the appropriate level for an arbitrarily chosen actuator loading. 
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Table 1 - g-LIMIT Parameters 


Parameter 

Symbol 

Value 




Flotor Mass 

m 

15.12 kg 

Flotor Moments of Inertia 


0.50 kg m 2 
0.62 kg m 2 
0.18 kg m 2 

Flotor Products of Inertia 


Ie-4 kg m 2 
-le-4 kg m 2 
-8e-4 kg m 2 

Umbilical Locations (F) 

(f v . 

[0.0 -0.12 -0.032] m 

(3 Umbilicals) 

f u 

[0.1 0.06 -0.032] m 

[-0.1 0.06 -0.032] m 

Actuator Current Vectors (S) 

(5) - 

[0.0 0.0 1.0] 

(6 Actuator Coils) 

Li 

[-1.0 0.0 0.0] 

[0.0 0.0 1.0] 

[0.5 0.866 0.0] 

[0.0 0.0 1.0] 

[0.5 -0.866 0.0] 

Actuator Magnet B-Field Vectors (F) 

(n B 

[0.0 1.0 0.0] 

(3 Actuator Magnets) 

— l 

[0.0 1.0 0.0] 
[0.866 -0.5 0.0] 

[0.866 -0.5 0.0] 
[-0.866 -0.5 0.0] 
[-0.866 -0.5 0.0] 

Actuator Constant 

( L , B,) 

1.0 N/Amp 


Table 2 - Actuator Current Pulse Values 



Scenario #1 

Scenario #2 


[Amp-sec] 

[Amp-sec] 

Actuator #1 Current Pulse 

0.0 

0.0 

Actuator #2 Current Pulse 

-0.231 

-0.10 

Actuator #3 Current Pulse 

0.0 

0.0 

Actuator #4 Current Pulse 

-0.139 

-0.10 

Actuator #5 Current Pulse 

0.0 

0.0 

Actuator #6 Current Pulse 

0.193 

-0.10 
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Figure 10. Flotor Energy for Scenario #1 with No Damping 
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Figure 1 1 . Flotor Energy for Scenario #2 with No Damping 
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Concluding Remarks 

This paper has presented the derivation of algebraic, state-space equations for the 
Glovebox Integrated Microgravity Isolation Technology (g-LIMIT). The states employed include 
payload relative translational position (jr n ) and velocity ( x b ), payload relative rotation (v (/ ) and 
rotation rate (x r ), and payload translational acceleration (*,.). Feedback of x (l corresponds to a 
change in effective umbilical translational stiffness, where the effective umbilical is assumed to be 
attached at the flotor center of mass. Similarly, feedback of x h , x d , or x e corresponds, 

respectively, to a change in effective umbilical translational damping, rotational stiffness, or 
rotational damping. Likewise, feedback of payload translational acceleration causes a change in 
effective payload mass. Thus, a cost functional which penalizes these states produces intuitive 
effects on system effective stiffness, damping, and inertia values. 

The acceleration states can be selected to pertain to any arbitrary point E on the flotor. 
This allows an optimal controller to be developed which penalizes directly the acceleration of any 
significant point of interest, such as the location of a crystal in a crystal-growth experiment. 

The equations have been put into state-space form so that the powerful controller-design 
methods of optimal control theory (e.g., H 2 synthesis, synthesis, // synthesis, mixed-// 
synthesis, and //analysis) can be used. References [7], [8], [13] and [14] detail the H 2 optimal 
controller design approach used for g-LIMIT, and Reference [15] describes the insights gained 
from a single-degree-of-freedom case study. 

The linearized state-space equations have been verified against a nonlinear truth model 
in a simulation study. The state responses from the model were shown to agree very well with the 
“truth” model for a set of parameters representative of the g-LIMIT configuration. Additionally, it 
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was shown that flotor kinetic and potential energy responses were as expected, demonstrating that 
the model conserves total energy for an arbitrarily chosen impulsive actuator loading. 
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